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1. Introduction 



o 

co 

C*~) ' Searching for a connection between classical and quantum integrable sys- 

tems is an old and new subject, occasionally leading to a breakthrough to- 
wards a new area of research. One of the landmarks in this sense is the quan- 
tum inverse scattering method, also known as the algebraic Bethe ansatz. 
Stemming from the classical inverse scattering method, the algebraic Bethe 
ansatz covers a wide class of integrable systems including solvable models 
of statistical mechanics on the basis of the Yang-Baxter equations. 1 More- 
over, remarkably, it was recognized later that a kind of classical integrable 
systems (discrete Hirota equations) show up in the so called nested Bethe 
ansatz. 2 

Recently a new connection was found by Foda and his group. 3-7 They 
observed that special solutions of the classical integrable hierarchies (the 
KP and 2D Toda hierarchies) are hidden in quantum (or statistical) inte- 
grable systems such as the 6- vertex model under the domain wall boundary 
condition (DWBC), 3 the finite XXZ chain of spin 1/2, 4,5 and some other 
quantum integrable systems. 6,7 Their results are based on a determinant 
formula of physical quantities, namely, the Izergin-Korepin formula for the 
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partition function of the 6-vertex model 8-10 and the Slavnov formula for the 
scalar product of Bcthe states in the XXZ spin chain. 11 ' 12 Those formulae 
contain a set of free variables, and the determinant in the formula is divided 
by the Vandermonde determinant of these variables. Foda et al. interpreted 
the quotient of the determinant by the Vandermonde determinants as a tau 
function of the KP (or 2D Toda) hierarchy expressed in the so called "Miwa 
variables" . 

In this paper, we review these results along with some additional infor- 
mation on this issue. We are particularly interested in the relevance of the 
2-component KP (2-KP) and 2D Toda hierarchies. Unfortunately, this re- 
search is still in an early stage, and we cannot definitely say which direction 
this research leads us to. A modest goal will be to understand the algebraic 
Bethe ansatz better in the perspective of classical integrable hierarchies. 

This paper is organized as follows. In Section 2, we start with a brief 
account of the notion of Schur functions that play a fundamental role in 
the theory of integrable hierarchies, and introduce the tau function of the 
KP, 2-KP and 2D Toda hierarchies as a function of both the usual time 
variables and the Miwa variables. Section 3 deals with the partition function 
of the 6-vertex model with DWBC. Following the procedure of Foda et 
al., we rewrite the Izergin-Korepin formula into an almost rational form 
and show that a main part of the partition function can be interpreted 
as a KP tau function. Actually, the partition function allows two different 
interpretations that correspond to two choices of the Miwa variables. We 
examine a unified interpretation of the partition function as a tau function 
of the 2-KP (or 2D Toda) hierarchy. In Section 4, we turn to the finite XXZ 
chain of spin 1/2, and present a similar interpretation to the scalar product 
of Bethe states (one of which depends on free variables) on the basis of 
the Slavnov formula. Section 5 is devoted to some other models including 
the phase model, 13 which is also studied by the group of Foda. 6 For those 
models, a determinant formula is known to hold for the scalar product of 
Bethe states both of which depend on free parameters. 14 We consider a 
special case related to enumeration of boxed plane partitions. 

2. Tau functions 
2.1. Schur functions 

Let us review the notion of Schur functions. We mostly follow the notations 
of Macdonald's book. 15 

For N variables x = (x\, . . . , Xn) and a partition A = (Ai, A 2 , . . . , \n) 
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(Ai > A 2 > • • • > Ajy > 0) of length (A) < N, the Schur function s\(x) can 
be defined by Weyl's character formula 

«M = A(x) ~ W 

where A(cc) is the Vandermonde determinant 

A(*)=det(z7^)£ =1 = TJ {xi-xj). 

l<i<j<N 

By one of the Jacobi-Trudi formulae, s\(x) can be expressed as a de- 
terminant of the form 

s A (x)=det(/i Ai _ i+j (x))5 =1 , (2) 
where h n (x), n = 0, 1, 2, . . ., are the completely symmetric functions 
K(x)= x kl x k2 ---x kn for n>l, fc (a:) = 1. 

l<fel<fc 2 <---<fcn<JV 

The complete symmetric functions h n (x) themselves can be identified with 
the Schur functions for partitions with a single part: 

h n (x) = S(„)(x), (n) := (n,0, . . . , 0). 

Another form of the Jacobi-Trudi formulae uses on the elementary sym- 
metric functions e n (x). Since we shall not use it in the following, its details 
are omitted here. 

The complete symmetric functions have the generating function 

oo N / N \ 

^2 K(x)z n = TJ(1 - x k zY x = cxp I -^log(l - x k z) I . 

n=0 fe=l V fe=l / 

Since log(l — x^z) has a Taylor expansion of the form 

oo n 

log(l - x k z) = - £ 

71=1 

this generating function can be rewritten as 

oo / oo \ 

^ & n (ss)3 n = exp I ^ i„z , (3) 

n=0 \n=l / 

where i n 's are defined as 

i„ = if>£. (4) 

n fc=i 
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In view of (2) and (3), one can redefine the complete symmetric functions 
and the Schur functions as functions h n [t] and s\[t] a of t — (ti,t 2 ,-.-), 
namely, 

oo / OO \ 

J2h n [t}z n = cxp (5) 
n=0 \n=l ) 

and 

SA [t] = det(/i Ai _ i+j [t])5 =1 , (6) 

where A is understood to be an infinite decreasing sequence A = (Ai, A2, . . .) 
with Ai = for all but a finite number of i's, and N is arbitrary integer 
greater than or equal to l(X) = max{i | \ ^ 0}. The right hand side of (6) 
is independent of TV, because the lower left block of the matrix therein for 
i > l(X) and j < l(X) is zero and the lower right block for i, j > l(X) is an 
upper triangular matrix with 1 on the diagonal line. 

This is a place where a connection with the KP hierarchy 17 shows up. 
Namely, the variables t = (t\,t2, ■ ■ •) are nothing but the "time variables" of 
the KP hierarchy, and the Schur functions s\[t] are special tau functions. As 
first pointed out by Miwa, 18 viewing the tau function as a function of the x 
variables leads to a discrete (or difference) analogue of the KP hierarchy. For 
this reason, the x variables are sometimes referred to as "Miwa variables" 
in the literature of integrable systems. 



2.2. Tau functions of KP hierarchy 

Let us use the notation r[t] for the tau function in the usual sense (namely, 
a function of t), and let t(x) denote the function obtained from r[t] by 
the change of variables (4). It is the latter that plays a central role in this 
paper. 

A general tau function of the KP hierarchy is a linear combination of 
the Schur functions 

T[t] = Y,C X 8 X [t], (7) 
A 

where the coefficients c\ are Pliicker coordinates of a point of an infinite 
dimensional Grassmann manifold (Sato Grassmannian). 19 Roughly speak- 
ing, the Sato Grassmannian consists of linear subspaces W ~ C N of a fixed 



a These convenient notations are borrowed from Zinn-Justin's paper. 
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linear space V ~ C z . We shall not pursue those fully general tau functions 
in the following. 

We are interested in a smaller (but yet infinite dimensional) class of tau 
functions such that c\ = for all partitions with l(X) > N. This corresponds 
to a submanifold Gr(N, oo) of the full Sato Grassmannian. The Pliicker 
coordinates c\ are labelled by partitions of the form A = (Ai, . . . , Ajv), and 
given by finite determinants as 

CA = dct(/ Mj )5 =1 , k:=\-i + N. (8) 

Note that the sequences Aj's and k's of non- negative integers are in one-to- 
one correspondence: 

oo > Ai > • • • X N > < — > oo > h > ■ ■ ■ > l N > 0. 

The N x oo matrix 

F = {fij)i=l,...,N, j=0,l,... 

of parameters represent a point of the Grassmann manifold Gr(A r , oo). 

By the Cauchy-Binet formula, the tau function t(x) in the a;-picture 
can be expressed as 

T[X) ~ , ^ , . A(x) _ A(x) ' { ' 

oc>h>--->l N >0 y ' y ' 

where /i(x)'s are the power series of the form 

1=0 

In particular, if there is a positive integer M such that 

fij =0 for i > M + N 

(in other words, fi(x)'s are polynomials of degree less than M + N), the 
Pliicker coordinate c\ vanishes for all Young diagrams not contained in the 
N x M rectangular Young diagram, namely, 

c A = for A % (M N ) := (M, . . . , M) 

S v ' 

JV 

The tau function r[t] thereby becomes a linear combination of a finite num- 
ber of Schur function, hence a polynomial in t. Geometrically, these solu- 
tions of the KP hierarchy sit on the finite dimensional Grassmann manifold 
Gr(N, N + M) of the Sato Grassmannian. 
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2.3. Tau functions of 2-KP hierarchy 

The tau function r[t, i] of the 2-component KP (2-KP) hierarchy is a func- 
tion of two sequences t — (ii, <2, • • •) and i = (ti, iz, . . .) of time variables, 
and can be expressed as 

r[i,t] =$> v s A [*K[*], (10) 

where ca^'s are Pliicker coordinates of a point of a 2-component analogue 
of the Sato Grassmannnian (which is, actually, isomorphic to the one- 
component version). 19 

The aforementioned class of tau functions of the KP hierarchy can be 
generalized to the 2-component case. Such tau functions correspond to 
points of the submanifold Gr(M + N, 2oo) of the 2-component Sato Grass- 
mannian. For those tau functions, the Pliicker coordinates vanish if 
> M or l(u) > N; the remaining Pliicker coordinates are given by finite 
determinants of a matrix with two rectangular blocks of size (M + N) x M 
and (M + N) x N as 

c X f, = det(/ Mj . | g iimk ), (11) 

where i is the row index ranging over i = 1, . . . , M + N and j, k are column 
indices in the two blocks ranging over j — 1, . . . , M and k = 1,...,N, 
respectively. Z/s and m^'s are related to the parts of A = (Xj)jL 1 and 
M = (Mi)jLi as 

= Aj - j + M, m k = Hk-k + N. 
By the change of variables from x and y to 



M N 

j=i *;=i 

the tau function r[t,t] is converted to the (a;, y)-picture r(x,y). Again 
by the Cauchy-Binet formula, r(x 7 y) turns out to be a quotient of two 
determinants as 

. , _ detCA(^) | gi (y fc )) 

1 ' yJ " A(aj)A(y) ' ( ' 

where the denominator is the determinant with the same block structure 
as (11), and fi(x) and gj(y) are power series of the form 



fi{x) = £ ^ lX ^ 9i{y) = ^2 guy 1 



1=0 
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2.4. Tau function of 2D Toda hierarchy 

The 2-KP hierarchy is closely related to the 2D Toda hierarchy. 20 The tau 
function r s [t,t] of the 2D Toda hierarchy depends on a discrete variable 
(lattice coordinate) s alongside the two scries of time variables t and t. 
For each value of s, T s [t,i] is a tau function of the 2-KP hierarchy, and 
these 2-KP tau functions are mutually connected by a kind of Backlund 
transformations. Consequently, r s [t,t] can be expressed as shown in (10) 
with the coefficients c s \^ depending on s. 

Actually, it is more natural to use s^[t]s„[— i] rather than s;Jt]s M [i] for 
the Schur function expansion of the Toda tau function. 21 (Note that s^[— i] 
can be rewritten as 

M-t] = (-i) H *vM, 

where t u denotes the transpose of u.) Expanded in these product of tau 
functions as 

the coefficients c s Xfi become Pliicker coordinates of an infinite dimensional 
flag manifold. Intuitively, they are minor determinants 

c sAm = det{g Xi -i+ s , H -0+3)^=1 (15) 

of an infinite matrix g = (gij)i,jez, though this definition requires justifi- 
cation. 21 In particular, if g is a diagonal matrix, the coefficients c s \^ arc 
also diagonal (namely, c s \^ oc 5\^) and the Schur function expansion (14) 
simplifies to the "diagonal" form 

00 

T s[t,i] = y^c s xs\[t]sx[-i], c s x = J^gxi-i+s- (16) 

If we reformulate the 2D Toda hierarchy on the semi-infinite lattice 
s > 0, the infinite determinants defining ca's are replaced by finite determi- 
nants, and T s [t,t] itself becomes a finite determinant. We shall encounter 
an example of such tau functions in the next section. 

3. 6-vertex model with DWBC 
3.1. Setup of model 

We consider the 6-vertex model on an TV x N square lattice with inho- 
mogeneity parameters u — (m, . . . , upf) and v — (vi, . . . , vn) assigned to 
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Fig. 1. Square lattice with DWBC 

the rows and columns. The boundary of the lattice is supplemented with 
extra edges pointing outwards, and the domain-wall boundary condition 
(DWBC) is imposed on these extra edges. Namely, the arrows on the extra 
edges on top and bottom of the boundary are pointing outwards, and those 
on the other extra edges are pointing inwards (see figure 1). 

The vertex at the intersection of the i-th row and the j-th column is 
given the following weight w^ determined by the configuration of arrows 
on the adjacent edges: 

— Wij = a{ui - Vj) 



w i: j = b(ui - Vj 



C(Ui - Vj) 



The weight functions a(u), b(u), c(u) arc defined as 

a(u) = sinh(u + 7), b(u) = sinhu, c(u) = sinli7, (17) 
where 7 is a parameter. Thus the partition function of this model is defined 
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as a function of the inhomogeneity parameters u and v: 

N 

Z N = Z N (u,v) = ^ H W ij- 

configuration i 'i= 1 



3.2. Izergin-Korepin formula for Zn 

According to the result of Korepin 8 and Izergin, 9 the partition function Zn 
has the determinant formula 

nfj=i sinh(u, - Vj + 7) smh( Ui - Vj) 

Zn — 



rii<i<j<iv sinh(u, - Uj) sinh(^ - v t ) 



- «le, ( — 7 -1" , (18) 

sinh(u, — Vj + 7) sinn(Uj — Vj) ' 



which one can rewrite as 

sinh 7 ^ X] 

Zn 



ni< l<7 <AT sinh( Ul - Uj) sinh(uj- - uj) 



Qfc=i sinh(Mj - ^ + 7) sinh( Mt - i> fc ) \ 
sinh(ui — + 7) sinh(ui — Vj) J ■ ■ 1 



and 

sinh w 7 



Z 



N 



ni<i<j<jv sinh(u, - uj) smh( Vj - Vi ) 

(N \ N 

IL=i sinh(Mfc - Vj + 7) sinh( Mfc - vj) \ 
sinh(wi — Vj + 7) sinh(ui — Vj) I 

If we introduce the new variables and parameters 3 

Xi \ e , yi . — e , q . e ^, 
we can rewrite these formulae as 

AT 

Zat = Cat J| (x^- 1 - yjq)(xi - 



1 j\ a. t q 1 - q 



X aTZTaT^ ** 7 * ~~T ) - (21) 



A(x)A(y) \(xiq 1 - y j q)(x l - yj ) 



N 
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Z N = C ^fj- f det ~ 7 fe) I , (22) 

A(x)A(y) [ (^q- 1 - yjq)(xi - yj) J ^ 1 ' 

and 

Z N = CN }f'- f dot ( U^fl - y)(** ZM Y , (23 ) 
A(x)A(y) \ (xiq^-yjqXxi-yj) . ' 

\ / i,j=i 

where Cat = Cn{u, v) is an exponential function of a linear combination of 
Ui's and Wj's. Apart from this simple factor, Zn thus reduces to a rational 
function of x = (xi, . . . , xn) and y = [y\, . . . , j/jv). 



3.3. KP and 2-KP tau functions hidden in Zn 

As pointed out by Foda et al., 3 two KP tau functions are hidden in these 
determinant formulae of Z N . Firstly, if y^s are considered to be constants, 
the cc-dependent part of (22) give the function 



V ' A(x) J3K ' (xq- 1 - y 3 q)(x - y 3 ) v ' 

This is a tau function of the KP hierarchy with respect to t n — i Y^k=i x k- 
Moreover, since /j(a;)'s are polynomials in x 1 this tau function is as a poly- 
nomial in t. In the same sense, if considered to be constants, the 
y-dcpendent part of (23) gives the function 

^M- gfr'sfK >M ,ifcy-^-», (25) 

A(y) (a;^ 1 - yq)(x l - y) 

which is a polynomial tau function of the KP hierarchy with respect to 

L n — n Z^fe=l ' 

Thus, apart from an irrelevant factor, Zn is a tau function of the KP 
hierarchy with respect to x and y separately. It will be natural to ask 
whether a tau function of the 2-KP hierarchy is hidden in Zn- 

A partial answer can be found in the work of Stroganov 22 and Okada. 23 
According to their results, if q = e 7 ™/ 3 , the partition function coincides, up 
to a simple factor, with a single Schur function of (x, y) as 

Zn = (simple factor)sA (x, y), (26) 

where A is the double staircase partition 



A = (N - 1, N - 1, N - 2, JV - 2, . . . , 1, 1) 
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of length 27V. By the way, for any partition A = (Ai, . . . , A27V) of length 
< 27V, the Weyl character formula for s\(x,y) reads 

det(^-^! I vt- l+2N ) 

SX{X ' V) = Mx^y) ' 

where the row index i ranges over i = 1, . . . , 27V and the column indices 
j, k in the two blocks over j, k = 1,...,7V. Multiplying this function by 
A(x,y)/A(x)A(y) gives 



A(x)A(y) AV A(x)A(y) 

which may be thought of as a 2-KP tau function of the form (13). 

Another answer, which is valid for arbitrary values of q, was found by 
Zinn- Justin (private communication). Let us note that (21) can be rewritten 

as 

JV 

Z N - CnIq- 1 - q) N ft (1 - q- 2 X iy f){\ - x^j 1 ) J[ (-y iVj ) 

i,j=l l<i<j<N 

( \ N 

x a / \ a / — tt det [ = =- | , (28) 

A(x)A( V -i) yil-q-HiyJ^l-Xiyj 1 )) V ; 

where 

y-^d/r 1 ,...,^ 1 ). 

The last part of thins expression, namely the quotient of the determinant 
by the Vandermonde determinants A(x)A(y _1 ), may be thought of as a 
2-KP tau function with respect to the time variables 

JV , TV 



n fc=l n fe=l 

This is a special case of the tau functions 



dct^^- 1 ))^ 
T(a; ^ )= A{x)A{y~ x ) ' (30) 

considered by Orlov and Shiota, 24 where 7i(z) is an arbitrary power series 
of the form 

oo 

h(z) = ^ Kz n , K ^ for n > 0. 

n=0 
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By the Cauchy-Binet formula, r(x,y) can be expanded as 

N 

r{x,y)= c x s x (x)s x (y- 1 ), c x := JJ h Xt - t+N - (31) 

A=(Ai,...,Aiv) t=l 

This is an analogue (for a semi-infinite lattice) of the Toda tau functions 
of the diagonal form (16). Note that the role of the lattice coordinate s is 
played by N. Since the number of the Miwa variables in (29) also depends on 
N, translation to the lanugage of the 2D Toda hierarchy is somewhat tricky, 
but this tricky situation is rather common in random matrix models. 24 
Thus, though not of the type shown in (13), the last part of (28) turns out 
to be a 2-KP tau function. 

Lastly, let us mention that Korepin and Zinn- Justin considered the par- 
tition function in the homogeneous limit as Uj, Vj — > 0. 25 In that limit, the 
partition function reduces, up to a simple factor, to a special tau function 
of the ID Toda equation, and can be treated as an analogue of random 
matrix models. 



4. Scalar product of states in finite XXZ spin chain 
4.1. L- and T -matrices for spin 1/2 chain 

We consider a finite XXZ spin chain of spin 1/2 and length N with in- 
homongeneity parameters £/, I = 1, . . . , N. To define local L- matrices, let 
us introduce the 2x2 matrix L(u) = (£ij(«))ij=i,2 of the 2x2 blocks 



Ln(u) = a(u) 



1 



b(u) 



1 



a 



L 12 {u) = c(u)a~ 



L 2 i(u) = c(u)a + , L 2 2(u) = b(u) 1 + ° + a(u) — 



where and cr 3 are the Pauli matrices 



1 






1 



1 

-1 



These 2x2 blocks are understood to act on the single spin space C 2 . The 
structure functions a(u), b(u), c(u) are the same as the weight functions (17) 
for the 6-vertex model, and built into the i?-matrix 



R{u - v) 



fa{u-v) \ 

b(u — v) c(u — v) 
c(u — v) b(u — v) 
V a(u-v)J 
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Let L( l \u — £;) = (L,fj(u — £))i,.j'=i,2 be the local i-matrix at the i-th 
site, 

L ij( u ~ £0 = • ' • ® 1 ® Lij{u - ® 1 <g> • • • , 
and define the T-matrix as 

n.)-(^»^»)- £ '"(.-ft)-^(.-£»). 

The matrix elements of these matrices, hence the trace of the T-matrix 

T(u) = TrT(u) = A(u) + D(u) 

as well, are operators on the full spin space V = C 2 . The //-matrices 

satisfy the local intertwining relations 

R(u - v)(L (l) (u) ® 7)(7 ® i (m) (w)) 

= (7 ® L^ m \v))(L^{u) ® J)i2(u - u), (32) 

where — u ), L(0( u ) eg) 7 and 7 <8> ^< m )(w) are understood to be 4 x 4 
matrices (of scalars or of spin operators on V) acting on the tensor product 
C 2 <8> C 2 of two copies of the auxiliary space C 2 . These local intertwining 
relations lead to the global intertwining relation 

R(u - v)(T(u) ® 7)(7 ® T(«)) = (7 ® T(«))(T(«) ® J)fl(u - v) (33) 

for the T-matrix. This is a compact expression of many bilinear relations 
among the matrix elements of T(u) and T(v), such as 

A(u)B(v) = f(u - v)B(v)A(u) - g(u - v)B(u)A(v), 

D(u)B(v) = f(u - v)B{v)D{u) - g(u - v)B(u)D(v), (34) 

C{u)B{v) = g{u - v)(A(u)D(v) - A(v)D(u)) 



and 



where 



[A(u),A(v)}=0, [B(u),B(v)] = 0, 
[C(u),C(v)]=0, [£>(«),£>(«)] =0, 



a(u) sinh(w + 7) c(u) sinh7 



(35) 



b(u) sinhu b(u) sinhw 

A consequence of those relations is the fact that T(u) and T(v) commute 
for any values of u, v: 

[T(u),T(v)] = 0. (36) 

The algebraic Bethe ansatz is a method for constructing simultaneous eigen- 
states (called "Bethe states") of T(u) for all values of u. 
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4.2. Algebraic Bethe ansatz 

Let us introduce the pseudo- vacuum 0} and its dual (0|: 

N N 



io)eV*, |o> = (g) ( I) gV. 



1=1 1=1 

They indeed satisfy the vacuum conditions 

(0\B(u) = 0, C(u)|0}=0, 

A(u)\0) = a(u)\0), D(u)\0) = S(u)\0), (37) 
(0L4(u) = a(u)<0|, (0|D(u)=J(«)(0|, 

where 

N N 

a(u) = ]^[sinh(u - & + 7), S(u) = J^sinh(u - 
;=i i=i 

For notational convenience, we introduce the reflection coefficients 

a(u) 

r(u) = Wj- 

Bethe states are generated from |0) by the action of B(vj)'s. Suppose 
that vj's satisfy the Bethe equations 

, , -r-r Sinh(«i — Vj — 7) . , 

^n S inhK-, jtt) =1 ' ,=1, -' n - (38) 

The state n"=i B( v i)\®) then becomes an eigenstate of T(u): 

n 

T(u)Y[B( Vl )\0) 

i=l 

(n n \ n 

«(«) n - m ) + n /( u - ^ n b k)i°>- (39) 
i=i i=i j i=i 

Let us remark that the operators A(u), B(u), C(u), D(u) are related to 
a row-to-row transfer matrix of the 6-vertex model on the square lattice. 
One can thereby derive the identities 10 

N 



(0\Y[C{u t ) = (0\Z N ( Ul ,...,u N , £ N ), 
i=l 

N 

Y[B{ Ui )\0) = Z N {ui,...,u N , !; N )\0), 



1=1 
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where Mi's are free (namely, not required to satisfy the Bethe equations) 
variables, and |0) and (0| denote the "anti-pseudo-vacuum" and its dual: 



N , Q 



<o| = (g)(oi)eV*, |0) = (g)r lev. 

i=i i=i ^ ' 

4.3. Slavnov formula for scalar product 

Let u = (u\, . . . , u n ) be free variables and v = (vi, . . . , v n ) satisfy the Bethe 
equations (38). The Slavnov formula 11 ' 12 for the scalar product 

n n 

s n (u,v) = (oin^K)n s ^)i°) 



i=l i=l 



reads 



rr?-i s(ui)5(vi) nr,-i smh(ui - « 7 - + 7) 

ni<i<j<„ smh(w 4 - U j) sinh(uj - Vj) 



where 



^7 [ ! _ r(u .) TT sinh(M t -Vk-j) 

13 sinh(ui - Vj + 7) sinh(u J - Vj) \ % ^ sinh(ui - v k + 7) 



One can rewrite this formula as 



where 



sinh(wi - Vk + 7) - a(t*») sinh(ui - w fe - 7) 



sinh(u, — Vj) 



4.4. JlP iau function hidden in S n (u,v) 

If we introduce the new variables and parameters 4 

Xi . — e % 1/^. e 1 , Zi . — e ^ , q . 6 ^ , 
the Slavnov formula can be converted to the almost rational form 
_ ^sinh^n^i^fa)^^^))"^! 
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where C n = C n (u, v) is an exponential function of a linear combination of 
ttj's and t>j's, and 

TV N 

]J(x - zi) Yliq^x - qy k ) - JJiq^x - qz { ) \\{qx - 
k^j 1=1 kjij 



fj(x) = ^— 

x - yj 

Thus, as pointed out by Foda et al., 4 a KP tau function of the form (9) 
is hidden in S(u,v). Moreover, since the Bethe equations (38) imply the 
equations 

N 



]J{yi - zi) Y[(q l Vi - qyk) 



1=1 



k^j 



N 



= Yl(l 1 Vi-QZi)Y[(qy i -q 1 y k ) (i = l,...,n) 

k^j 



1 = 1 



for yj's, the numerator of fj(x) can be factored out by the denominator 
x — yj. Thus fj(x)'s turn out to be polynomials in x. 



5. Scalar product of states in models at q — 

A class of solvable models, such as the phase model 13 and and the totally 
asymmetric simple exclusion process (TASEP) model, 26 can be formulated 
by a set of 2 x 2 L-matrices (u), I = 1, 2, . . . , N, and an i?-matrix of the 
form 



R(u - v) 



(f{u-v) \ 

1 g(u-v) 

g{u-v) 

\ f(u-v)J 



where 



f(u-v) 



l(u - v) 



uv 



This i?-matrix is obtained as a "crystal" (namely, q — > 0) limit of the R- 
matrix of the 6- vertex model and the XXZ spin chain. Unlike the XXZ chain 
of spin 1/2, the L-matrices are not given by 2 x 2 blocks of the i?-matrix 
and take a mo del- dependent form. We define the T-matrix as 

™=(cS:!^)= i(i, w- i< " , w 
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and consider the scalar product 

n n 

S n (u,v) = 

»=i »=i 

Remarkably, even if both u = (ui ,...,«„) and v = (v\, . . . ,v n ) are free 
variables, the scalar product for those models has a determinant formula 14 
of the form 

S n {u,v) = (simple factor) [j 2 M '_" j 2 ^ 2 det(^-)^ =1 , (43) 

i<i<j<„ u * ~ U J U J ~ ^ 

where 

_ atu^t^K^M)"" 1 - ^(uQafa-)^/^)"- 1 
(«? - 

a(«) and r5(u) being determined by the action of A{u) and D(u) on the 
pseudo-vacuum. If a(u) and S(u) are given explicitly, we will be able to 
obtain a KP tau function (and hopefully a 2-KP tau function as well). 

Such an interpretation can be found most clearly in the cases of the 
phase model 27-29 and the 4-vertex model, 30 both of which are related to 
enumeration of boxed plane partitions. The scalar product of Bcthc states 
in these cases becomes, up to a simple factor, a sum of products of two 
Schur functions: 

S(u, v) = (simple factor) s\{u\, . . . ,u 2 n )s\{v^ 2 , . . . ,v~ 2 ). (44) 

AC(JV") 

By the same reasoning as the interpretation of (31), one can see that this 
sum is a tau function of the 2-KP hierarchy (or, rather, the 2D Toda hier- 
archy as Zuparic argued 6 ) with respect to the time variables 

i=l i=l 

Actually, this case admits yet another interpretation. In the course of 
deriving (44), the scalar product is shown to be a generating function for 
counting plane partitions in a box of size n x n x N. It is more or less well 
known 16 that the Schur function S( N ^ gives such a generating function. 
Thus, in terms of this Schur function, the scalar product can be expressed 
as 

S(u,v) = (simple factor)s(7v«) (u 2 , . . . ,u 2 N ,v\, . . . ,v%). (45) 
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This is similar to the formula (26) of the partition function of the 6- vertex 
model for q = e 71 '/ 3 ; one can thereby find an interpretation as a KP or 2-KP 
tau function with respect to the time variables 

1 n 1 n 

n *r~i n 
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